ON TWO-DISTILLABLE WERNER STATES 



DRAGOMIR Z. DOKOVIC 



Abstract. We consider bipartite mixed states p in a d Cgi d quantum system 
with d > 3. We say that p is PPT if its partial transpose 1 ® T{p) is positive 
semidefinite, and otherwise p is NPT. The well-known Werner states are par- 
titioned into three types: a) the separable states (same as the PPT states), 
b) the 1-distillable states (necessarily NPT), and c) the NPT states which are 
not 1-distillable. We give several different formulations and provide further 
evidence for validity of the conjecture that the type c) Werner states are not 
2-distillable. 



1. Introduction 

Let Ti, — Ha'^'Hb be the Hilbert space for the quantum system consisting of two 
parties, A and B (Ahce and Bob). We assume that the Hilbert spaces Ha and Hb 
have the same finite dimension, which we denote by d. A product state is a tensor 
product Pa® Pb of the states pA and of the first and second party, respectively. 
A bipartite state p is separable if it can be written as a convex linear combination 
of product states. We say that a bipartite state is entangled if it is not separable. 
We say that p is PPT if its partial transpose a — 1 ®)T{p), computed in some fixed 
orthonormal (o.n.) basis of H_b, is a positive semidefinite operator. Otherwise a 
has a negative eigenvalue, and we say that p is NPT. 

It is more complicated to give the definition of distillability for bipartite states 
p. For that purpose we have to consider multiple copies of p. For k copies the 
density matrix is the fcth tensor power p*^*^ which acts on the Hilbert space H®'^. 
We can identify n'^'' with the tensor product of the Hilbert spaces nf" and T^f". 
In this way we can view p®'^ as a bipartite state. Thus any vector 1-0) G 'H'^*' has 
its Schmidt decomposition and a well-defined Schmidt rank. 

The definition of distillability given below is not the original one but it is the 
only one that we are going to use. Replacing the original definition with this one 
was nontrivial, see [5]. 

Definition 1.1. For a bipartite state p acting on % and an integer k > 1, we say 
that p is fc-distillable if there exists a (non-normalized) pure state j?/') € 'H®'^ of 
Schmidt rank at most two such that 

(1.1) (V'la^'^IV) <0, 

where a — 1® T{p) is the partial transpose of p. We say that p is distillable if it is 
k-distillable for some integer k>l. 
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If a bipartite state p is separable then it is PPT, i.e., a is positive semidefinite, 
and consequently p is not distillable. For the same reason, the entangled bipartate 
PPT states are not distillable. Equivalently, every distillable bipartite state is 
necessarily NPT. It is not known whether the converse holds, i.e., whether every 
bipartite NPT state is distillable. However it is widely believed that the converse 
is false. In particular the following conjecture has been raised [2]. 

Conjecture 1.2. There exist bipartite NPT states which are not distillable. 

It is known [10 that for each integer fc > 1 there exist examples of bipartite 
states which are distillable but not fc-distillable. 

We fix an o.n. basis i = 1, . . . , d of Ha, and an o.n. basis of "Hb for which 
we use the same notation. The context will make clear which basis is used. After 
fixing these bases, we can define the flip operator _F : "H — >■ "H by 

The (non-normalized) Werner states on % can be parametrized as 

(1.2) pw{t) = l-tF, -l<t<l. 

Let IV'max) € 'H be the maximally entangled (pure) state given by 

Its density matrix is the projector 

Since dP is the partial transpose of F, the partial transpose of pw{t) is 

awit) = 1 - tdP. 
The following facts about the Werner states are well-known. 

Proposition 1.3. The Werner states pw{t) are: 

(a) separable if and only if ^1 < t < 1/d; 

(b) NPT if and only ifl/d<t<l; 

(a) l-distillable if and only if 1/2 < t < 1. 

The importance of Werner states for the distillability problem for bipartite states 
was first established in [4 . From now on we assume that d > 3. 

Proposition 1.4. Conjecture \1.2\ is equivalent to the assertion that some NPT 
Werner states pw{t) are not distillable. (Such t necessarily belong to the interval 
(lMl/2].j 

In fact the following stronger conjecture is believed to be true [UEllS]. 
Conjecture 1.5. None of the Werner states pw{t), t E (l/d, 1/2], is distillable. 

In this paper we will consider a very weak version of it. 
Conjecture 1.6. None of the Werner states pw{t), t E {1/d, 1/2], is 2 -distillable. 

For the fc-distillability problem the following fact ^ is important. 
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Proposition 1.7. If pw{to) is not k-distillable for some to, then none of the states 
pw(t), t <to, is k-distillable. 

In view of this proposition, it suffices to prove Conjecture 11.61 for t = 1/2 only. 
Extensive numerical evidence for the validity of this conjecture in the case d = 3 
is presented in [2l |3l [7] and [9j . In the latter it is also claimed that their numerical 
proof is rigorous. The case d = 4 was analyzed recently in [8 but it remains open. 

The paper is organized as follows. In Section [2] we construct a hermitian bi- 
quadratic form $ and show that Conjecture 11.61 is equivalent to $ being positive 
semidefinite, $ > 0. The form $ depends on Ad arbitrary vectors Xi,yi S "Ha and 
e Hb, i = I,.. ■ ,d. 

In Section [3] we express <& as a function of four matrices X, Y, U, V of order d, 
where X = [xi, . . . , Xd], etc. The formula that we obtain for $ is very attractive 
and looks simple, but proving that $ > remains a challenge. 

In Section |4] we compute the matrix H oi ^ when the latter is viewed as a 
hermitian quadratic form in the 2d^ complex entries of U and V. The entries of X 
and Y play the role of parameters. As expected. Conjecture [T|6] is equivalent to H 
being a positive semidefinite matrix. 

In Section [5] we prove that H{X,Y) > when X and Y are diagonal matrices. 
We point out that H{X, Y) is not diagonal even when both X and Y are. Since this 
is done for arbitrary d > 3, and the proof is quite nontrivial, this is an important 
piece of evidence for the validity of Conjecture 11.61 

In Section [6] we propose much simpler conjecture. Namely we replace the con- 
dition H{X,Y) > by just det H{X,Y) ^ with a mild condition on the pair 
{X,Y). The condition says that the pair {X,Y) is generic which means that the 
matrices X and Y are linearly independent and that some linear combination of 
them is nonsingular. We prove that this new conjecture implies Coniecture ll.61 but 
we failed to show that they are equivalent. 

2. The hermitian biquadratic form $ 
Since we are going to use only one Werner state, the one for t — 1/2, we set 
Pw = Pw(l/2) - 1 - F/2, <7w = ow{\l2) = 1 - dP/2. 
Conjecture II .61 is equivalent to the claim that the inequality 



is valid for all S H®'^ of Schmidt rank < 2. Such 1-0) can be written as 

Wl = \i>\) + 1-02), where 



Note that \x), \y) G V.a^T~1-a while \u), \v) e 'Hb®'-Hb- We point out that we do not 
require + 102) to be the Schmidt decomposition of \'ip)^ i.e., we do not require 
that {x\y) = {u\v) = 0. The reason for this is to allow the vectors |a;), \y), \u), \v) to 
be completely arbitrary. 

We can rewrite j^i) and 102) as 



(2.1) 



(V'|a|.2|0)>O 



|0l) = |x) (g) |u), \^^) = \y)®\v). 
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The vectors \xi) and \yi) live in Alice's second copy of Ha, while \ui) and \vi) live in 
Bob's second copy of "Kb- The summation is taken over all i and j in {1, 2, . . . , d}. 
Consequently, we can view the LHS of (|2.ip as a function of Ad vectors Xi, yj, Ur,Vs: 

. . . ,Xd,yi, . ■ . ,yd,ui, . . . ,Ud,vi, . . . ,Vd) = {ip\a^\ip). 

As 

a®2 = 1 - 1(1 (g) dP + dP (S) 1) + ^dP (g) dP, 

we have 

1 . X 1 

where 



$1 = (^1^), 

$2 = {Ml(E>dP\'iP), 

$3 = (^|dP^l|^), 

$4 = (V'ldP^dPIV'). 

After the substitution = + |'02), each of the $fe breaks up into four pieces. 
For instance, we have 

$2 = {i,j,Xi,Uj\l®dP\r,S,Xr,Us) 
i,j.r,s 

+ X! {'^d,Xi,Uj\l^ dP\r,s,yr,Vs) 

i.j,r,s 

We have computed each of the resulting 16 pieces. For instance the second piece, 
say E, in the above formula for $2 is computed as follows. We first observe that 
(i, J, Xi, Uj\l (El dP\r, s, yr, Vg) — Oii r^i or s^j. Thus we have 

E = ^{xi,Uj\dP\yi,Vj) 

i.j \ r s J 

= ^(:r,|u*)(y,|z;*)*. 
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The final formulae are: 





$2 



+ ^(y.|t,;)(a:.K)*+^|(2/.|«;)r, 




$3 



^ (xi I ) (-u, I ) + ^ (x, I y, ){u,\vj) 





$4 




These formulae show that each <i>fe , viewed as a function of the components of 
the Xi and y^, is a hermitian quadratic form. And the same is true when we view 
them as functions of the components of the Ui and Vj . Hence we shall refer to the 
(and $) as hermitian biquadratic forms. We note that $ remains a hermitian 
biquadratic form when we replace the partition of the vector variables Xi, yj\ur, Vs 
with the new partition Xi,Vs\yj,Ur. (This is a consequence of the fact that the 
two-pair Werner state is invariant under the flip operator 1 (g) F of the second 
pair.) 

The following proposition follows immediately from Eq. (j2.ip and the definition 
of the form $. 

Proposition 2.1. Conjecture \1.6\ is equivalent to the assertion that the form $ is 
positive semidefinite. 



It is quite natural to introduce the d x d matrix X whose successive columns 
are the vectors xi, . . . , Xd- Similarly, we define the matrices Y, U and V . Let 
denote the space of complex matrices of order d. Define the inner product on Md 
by {A\B) — tr (A^B). For the corresponding norm we have ||j4|p = tr (A^A). 
The tensor product of matrices A = [a^] and B is defined as the block-matrix 
A(^B = [aijB]. 

Now the formulae for $ can be rewritten in terms of the matrices X, Y, U and 
V. We obtain that 



^i{X,Y,U,V) = \\Xf\\Uf + \\Yf\\Vf + 2^{tr{X^Y)-tr{U^V)), 

<!>2{X,Y,U,V) = \\X^U + Y^Vf, 

<!>3{X,Y,U,V) = tr {X^X*U^U + X^Y*V''U + Y^X*U''V + Y^Y*V^V) , 

<S>4{X,Y,U,V) = \tr {X^U + Y^V)\'^, 



3. $ AS A FUNCTION OF FOUR MATRICES 
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where 3? stands for "the real part of" . 

The first expression can be further simplified by using the standard Frobenius 
norm on the tensor product of matrices 

Y,U,V) ^\\X (E)U + Y (E) Vf. 

The third expression also simplifies to 

^3{X,Y,U,V) = \\UX^ + VY'^f. 

Consequently, we have 

^{X,Y,U,V) = \\X (S)U + Y (S)V\\^ ~^\\X'^U + Y^V\\^ 

^l\\UX^ + VY^W' + i|tr {X^U + Y^V)\^. 

One benefit of these formulae is that they make the verification of the next 
proposition trivial. 

Proposition 3.1. The hermitian biquadratic forms $fc (and<^) are invariant under 
the following two types of transformations: 

(a) X ^ AX, Y ^ AY, U A*U, V A*V; 

(b) X^XB, Y^YB, U^UB*, V^VB*, 

where A and B are arbitrary unitary matrices. Thus we have an action of the group 
U(d) X \J{d) on the matrix quadruples [X, Y, U, V) which preserves these forms. 

We can now reformulate Proposition 12. II as follows. 

Proposition 3.2. Conjecture \1.6\ is equivalent to 

\\X®U + Y® V\\^ - ^WX'^U + Y^Vf - ^\\UX'^ + VY'^lf 

+ i|tr {X^U + Y'^V)\^ > 0, \fX,Y,U,V e Md. 

Let us remark that in view of the previous proposition we can assume, without 
any loss of generality, that the matrix X is diagonal with real nonnegative diagonal 
entries. 



4. The matrix H of the form $ 



We shall consider the entries of X and Y as parameters and those of U and 
V as complex variables. Then $ (and each becomes a family of hermitian 
quadratic forms depending on the mentioned parameters. Since we have derived 
explicit formulae for $, it is a routine matter to compute the matrix H = H{X, Y) 
of $. For that purpose we used the formulae given in Section [2j We shall only 
give a very simple recipe for the construction of H. As $ depends on 2d^ complex 
variables, its matrix is a hermitian matrix of order 2d^. We shall write for the 
identity matrix of order m. 

For any complex matrix Z let Z denote the column vector obtained by writing the 
columns of Z one below the other starting with the first column, then the second, 
etc. Now we can express the relationship between the form $ and its matrix H by 
the formula 



(4.1) 



^X,Y,U,V) = 



W ft 



H{X,Y) 



U 
V 
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Similar formulae are valid for the forms $fe and their matrices Hk- 
The matrices Hk are given by the following simple formulae: 

tr(xty) 
tr(rtx) ||y||2 

X^X X^Y 
Y^X Y^Y 

X*X'^ Id < 
-Y*XT 

X*X^ ^*yT ■ 

Y*xT y*YT 



Id' 
Id' 



.X*Y^ 
Y*Y^ 



Those for Hi and H^ are obvious. We omit the tedious but straightforward verifi- 
cation of the formulae for H2 and H^. 
For H we have the formula 



(4.2) 



H{X,Y) = Hi-]^{H2 



H. 



In view of Proposition [2Tl we can restate Conjecture 1 1 . 61 in the following equiv- 
alent form. 

Conjecture 4.1. The matrix H{X,Y) is positive semidefinite for all X^Y ^ M^. 

Recall that d > 3 by the assumption made earlier, but this conjecture is also 
valid for d = 1 and d = 2. However in these two cases the determinant of H{X, Y) 
is identically 0. 

If A € U((i) and we replace X and Y with AX and AY, respectively, then the 
Hk undergo the transformation Z {l2d ® A*)Z(l2d ® A'^). In fact Hi and H2 
remain fixed under this transformation. 

Similarly, if _B e V{d) and we replace X and Y with XB and YB, respectively, 
then the H^ undergo the transformation Z — > {I2 ^ B'' Id)Z{l2 i? (8) Id)- This 
time Hi and H^ remain fixed. 

Hence we have the following formula. 

Proposition 4.2. For A, B E V{d) we have 

(4.3) H{AXB, AYB) = {h ® B'^ ® A*)H{X, Y){h ®B® A^). 

This formula is also valid for the matrices Hk ■ 

Next we show that H{X,Y) satisfies yet another identity. Let 

a /3 



A 



7 



(A^)- 



a 
7' 



/3' 
5' 



(4.4) 
and 
(4.5) 

By using 



we deduce that 

\ax + Py) (g) \a'u + P'v) + + Sy) ® \^'u - 



e GL2(C) 



a 


7 




'a' P' ' 




■ 1 


" 


. P 


5 




_ i 5' 







1 



-5'v) = \x) ® \u) + \y) ® \v) 
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and consequently 

$(aX + PY, -fX + SY, a'U + P'V, 7'C/ + 5'V) = F, U, V). 
By invoking Eq. ()4.ip and by using the formula 




we deduce that 



(4.6) H{aX + I3Y,^X + 5Y) = (A* (g) /<j2)H(X, y)(A^ ® /<j2). 



5. The diagonal case 



We say that the pair of complex matrices (X, Y) in Md is generic if X and Y 
are linearly independent and some linear combination of X and Y is nonsingular. 
In this section we analyze the structure oi H = H{X,Y) for a generic pair {X,Y) 
of diagonal matrices and prove Conjecture 14.11 in that special case. We denote the 
diagonal entries of X and 1" by Ai, . . . , and fj,i, . . . , fid respectively. 

Theorem 5.1. If {X,Y) is a generic pair of diagonal matrices, then H{X,Y) is 
positive definite. 

Proof. After replacing H with IITJII"^ where 11 is a sutable permutation matrix, H 
becomes direct sum of — d blocks of order 2 and an additional block of order 2d. 
It suffices to show that each of these blocks is positive definite. 

The blocks of order 2 are indexed by the integers p = (i — l)d + j, where 
i, j € {1, . . . , d} and j ^ i. For such index p, the corresponding block of order 2 is 
the principal submatrix H{p) of the original matrix H corresponding to indexes p 
and p -\- . Explicitly we have 



where — 1 for k ^ i,j and Ci — Cj — 1/2. Each matrix on the RHS is positive 
semidefinite. If H{p) is singular, then all of these matrices must be singular and 
must have the same kernel. This contradicts the linear independence of X and Y. 
Hence H{p) must be positive definite. 

It remains to consider the remaining block B of size 2d, i.e., the principal subma- 
trix of H corresponding to the indexes {i~ l)d + i and {i — l)d + i + d'^ for 1 < i < d. 
We have B = Bi — {B2 + B^)/2 + B^^/A where Bk denotes the corresponding prin- 
cipal submatrix of H^. Let us first consider the matrix B' — Bi ~ {B2 + Bj)/2. 
After a suitable simultaneous permutation of rows and columns, B' breaks up into 
the direct sum of d blocks G{i) of order 2, where i £ {1, . . . ,d]. Explicitly we have 



Each G{i) is positive semidefinite of rank 1 or 2. Thus in the decomposition 
B — B' + B4 both summands, B' and B4, are positive semidefinite. Assume 
that, say G(l) is singular. Then all d — 1 matrices in the above formula for G{1) 
must be singular and have the same kernel, i.e., the vectors (/i2, Ms, ■ • ■ , /id) and 
(A2, A3, . . . , Xd) are linearly dependent. Since X and Y are linearly independent, 
all the other G{i) must be positive definite. We deduce that the kernel of B' is 



d 



>'ktil lAifcP 
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1-dimensional and is spanned by the column vector having all components except 
the first which is ~/i2 and (d + l)st which is A2. Since this vector is not killed by 
B4, we conclude that B is positive definite. □ 

Corollary 5.2. Conjecture \4.l\ is valid when X and Y are diagonal matrices. 

Proof. This follows from the theorem because any pair of diagonal matrices can be 
approximated by a generic pair of diagonal matrices. □ 



6. A SIMPLIFIED CONJECTURE 

We have already observed that in proving that ^{X,Y,U,V), or equivalently 
its matrix H{X,Y), is positive semidefinite we may assume that X is diagonal 
with real nonnegative diagonal entries. In fact we may clearly assume that these 
diagonal entries are positive. As the size of H is rapidly increasing with d (it is 18 
for 0? = 3) and it depends on many parameters, the task of proving that H is positive 
semidefinite seems formidable. The standard criterion would require checking that 
all principal minors of H are positive semidefinite. For positive definiteness of H 
the leading principal minors would suffice. 

Let us partition H into four square blocks of order d^. Then it is not hard to 
show that the two diagonal blocks are positive semidefinite. In the case d ~ 3 
we were able to show that also the leading principal minor of H of order 10 is 
a positive semidefinite polynomial. This was accomplished by finding an explicit 
expression for this minor as a sum of squares of real polynomials. Although some 
minor further simplifications are possible, it is clear that such a brute force method 
will not provide a proof of the conjecture in general. Even in the case d = 3 we 
failed to finish the job. 

We propose now a new simpler conjecture which is stronger than Coniecture l4.1l 
We introduce the real valued polynomial D{X,Y) — det H{X,Y). By taking the 
determinants in Eq. (|4.3|) we deduce that D{X,Y) is invariant under the action of 
the group lJ{d) x lJ{d), i.e., 

(6.1) D{AXB, AYB) = D{X, F), VA, B e U(d). 
From Eq. (14. 6p we deduce that 

(6.2) D{aX + PY, + 6Y) = \aS ~ Pj\'^'^^D{X, Y) 

is valid when A is invertible. Since both sides are polynomials, this identity must 
be valid for arbitrary A. 

Note that D{X,0) = for all matrices X. More generally, we claim that 
D{X,Y) = if A and Y arc linearly dependent. Indeed, it suffices to choose a 
matrix A as in (|4.4p such that 7A + SY = and apply Eq. (|6.2p . The converse of 
this claim is false but we conjecture that it is true in a weaker form. 

Conjecture 6.1. D(X,Y) ^ for generic {X,Y) and d > 3. 

Proposition lS.ll shows that this conjecture is true when the matrices X and Y are 
diagonal. As this conjecture deals with only one polynomial and has no positivity 
conditions whatsoever, it should be much easier to prove (or disprove). 

Proposition 6.2. Conjecture \4.1\ is a consequence of Conjecture \6.1\ 
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Proof. Let Xi and Yi be any matrices in Md- We have to show that H{Xi,Yi) 
is positive semidefinite. Clearly it suffices to prove this when the pair {Xi,Yi) is 
generic. Let {XqjYq) be a generic pair of diagonal matrices. Then H{Xo,Yq) is 
positive definite by Theorem 15. II Consequently, D(Xo,Yq) > and all eigenvalues 
of H{Xq, Yq) are positive. We can join the pairs {Xq, Yq) and {Xi, Yi) by a continu- 
ous path {Xt, Yt), < t < 1, such that {Xt, Yt) is generic for each t. By Conjecture 
16.11 D{Xt,Yt) for all t. Hence H{Xt,Yt) has no zero eigenvalues. Since the 
eigenvalues of H{Xt,Yt) are continuous functions of t and they are all positive for 
t — 0, they must all remain positive for all values of t. In particular this is true for 
t = I. We thus conclude that H{Xi,Yi) is positive definite. □ 

We do not know whether Conjectures 14.11 and 16.11 are equivalent. 
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